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A new reconstruction method for Wigner function is reported for quantum tomography based on compressed 
sensing. By analogy with computed tomography, Wigner functions for some quantum states can be recon- 
structed with less measurements utilizing this compressed sensing based method. 
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Wigner function(WF), a quasi-probability distribution in phase space, was first introduced to describe quantum state in quan- 
tum mechanics by E. P. Wigner[r]. And later, it was extended to classical optics and signal processing. Since its birth, a 
great number of applications have been conducted in different fields. As for the original quantum case, like a continuous one- 
dimensional quantum system, the Wigner function is defined as 
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where x and ^ are positions and p is momentum, p denotes the density operator, and h is set to 1 for simplicity. 

One of the advantage for the Wigner formalism is that a tomographic scheme which could be used to reconstruct the WF. 
On the other hand, in medical tomography, like X-ray computed tomography. X-ray photons transmitted through the patient 
,— — , along projection lines from one side, and detectors measured the number of transmitted photons on the other side, then from 
the distribution of X-ray photons in different incident angles, linear attenuation coefficients in the slice being imaged are re- 
Q-i covered. Likewise, in phase space tomography, the WF tomographic reconstruction is more or less the same, based on a set of 
intensity measurements! 2]. This is because both of the two tomographic schemes share the same mathematical basis — Radon 
transform(Eqn.(2)) and its inverse(Eqn.(3)). 
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' Where P denotes the Cauchy's principle value, pr{x^,4>) denotes the measured distribution, in f{x,y) is the quantity in to- 
mography to be recovered. In tomography, objective function is reconstructed using back-projection algorithm for the inverse 
Radon transformation(Eqn.(3)). However, such a algorithm has some drawbacks, for example, many measurements are needed 
^— I to achieve a relatively good recovery. To reduce the dose of radiation on the patient in CT or make a more efficient estimate 
of quantum state in phase space tomography, the amount of measurements should be as small as possible but guarantee a good 
enough recovery. Nevertheless, according to the well-known Nyquist-Shannon sampling theorem^, it seems impossible to 
depress the measurements without giving a negative impact on the resolution. 

Recently, with the arise of a novel theory called compressed sensing (compressive sampling or CS) 101, people find a funda- 
mentally new approach to acquire data. It demonstrates that if the images or signals are sparse, we can recover them from what 
was believed to be highly incomplete measurements. In fact, it is just the application to medical tomography that started the 
seminal paper in CS. For example, in CT, X-ray images of human body are sparse in the wavelet and gradient representation, 
then compressed sensing is applied to reduce the irradiation and obtain clearer medical images fs*, '6'l. And then an immediate 
question would be is it possible to reconstruct WF in phase space tomography via compressed sensing? In fact, there exists 
infinite family of quasi-probability distribution in phase, including Husimi-Kando Q-Function, Glauber-Sudarshan P-Function 
and so on, but WF is chose for it is directly related to the quadrature histograms and can be easily measured through experiments. 

In this brief report, we explore the similarity between phase space tomography and computed tomography, and report a new 
discrete WF reconstruction method based on compressed sensing. We mention that some compressed sensing protocols have 
been proposed to reconstruct density matrix though quantum tomography [Til. 

Suppose that a signal x E is sparse in some basis x — 'i'x', which means x' has k{k <C N) nonzero entries. Then a 
measurement matrix <i> E (£mxN jg n^Qd to sense x and obtain a measurement vector y E C™. The central idea of compressed 
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sensing is that it is possible to reconstruct sparse signals of scientific interest accurately and sometimes exactly by a number of 
incoherent sampling which is far smaller than N. In general, such a recovery of x by performing ^o-minimization is a NP-hard 
problem. An alternative procedure called i!i -minimization is usually used which says 



min ||a;'||f^ subject to y ^ ^^ifx' ^ Tx' , (4) 

This is a convex optimization problem, and many numerical algorithms apply for its solution. 

Before showing the reconstruction method, we note that there are two main differences between conventional medical to- 
mography and the quantum tomography. First, WF is a pseudo probability distribution, may have negative values, while the 
attenuation coefficients in CT images are always nonnegative. Second, for a given angle, a complete probability density for the 
quadrature operator is measured due to the quantum mechanical property, but in medical tomography, detectors' number can be 
adjusted. On the other hand, similar to the medical case, at the heart of a possible CS application to phase space tomography is 
to find a suitable sparse representation basis. In real computed tomography, images are generally extended distributions which 
violate the prerequisite of ^i-based algorithms. But other sparseness property (gradient sparse) and reconstruction algorithm(TV- 
algorithm)|5] are found to reconstruct certain tomographic images which are relatively constant over extended areas. While due 
to its configuration's flexibility , it is hard even impossible to find a unitive sparse basis for WF in general. However, for some 
special applications, WF might share the same sparse basis. And in contrary to the medical images, some quantum states have 
finite values on a small region in phase space with zero or nearly zero values elsewhere, so have been sparse in this pseudo- 
pixel representation. Other situations may be completely different, for WF is not sparse in phase space representation. A basis 
transform like Fourier transform, wavelet transform may bring WF to a sparse representation. 

Although it is also possible to reconstruct continuous Wigner function, but discretization must be taken before it can be tackled 
in practise. So in this report we adopt a discrete quantum tomography model. But we note that it should be straightforward to 
transport our idea to an continuous one. And indeed such an attempt has been taken in information sciencel6]. The quantum 
state tomographic model for discrete Wigner functions first proposed by U. Leonhardt lfTlTl . We note that continuous and discrete 
quantum tomography are the same in principle. In finite quantum systems, position and momentum take discrete values in Z^. 
In Leonhardt's discrete quantum tomography picture, a WF in a (i-dimensional(odd) Hilbert space was proposed(Eqn.(6)) and 
phase was quantized to take place of momentum in continuous variable quantum systems. 
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Here, m, n denotes the spins or particle numbers, fi, v denotes the quantized phases. Similar to the continuous case, the observ- 
able distribution is the overlap of measured state's WF and the observable state's WF in phase space. 

pr{f^o,T) = ^ W(m, fi)S{n - Mo + 2mr; d) 

pr(mo, r) — W{m, ii)d{m — mo; d)6{m — mo; a). (6) 

Here, mo, Ho G Zrf,pr(/io, r) denotes the probability distribution, W{m, ^) denotes the WF to be recovered, r £ Z^, m, niQ 
represent the spins or particle numbers, ^, no represent the quantum phases and S{k; a) is the modular Kronecker symbol(equals 
1 when fc=0 (mod a), and zero otherwise) 

In our simulation, a WF for finite dimensional coherent state ||a|e"^)s (s=d-l) lfl2ll is used as an example to illustrate the CS 
based reconstruction algorithm. 



k=M+l 

exp[i(2fc - A/ - s - l)(27rM/d-(/) + 7r/2)] 



k=0 



[/c!(M-/c-s - l)!]i/2 

-M){2TTn/d- 
[fc!(Af-fc)!]i/2 



Gik 



^ cxp[»(2fc - M){2nfi/d -cj) + tt/2)] ^ 

2^ fhuw^mm 



3 



liii 



m 



i 



i 



(a) 



(b) 



FIG. 1: Darker regions correspond to higher values, the regions between dark and white are dominated by zero or nearly zero values 
.(a)Accurate WF for finite dimensional coherent state(Eqn.(8), d — 19, |a — 1.472)). (b) Recovered WF via compressed sensing, 285 
rows have been taken from total measurement matrix as the sensing matrix. 
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with 7] = 0, 1, M = 2nmod((i), Y\ts+i{x) is the Hermite polynomial (He„(a;) = 2 "/^i7„(a;/-\/2)) and is the fcth root of 
He,+i(x) =0. 

It should be pointed out that different from classical signal processing (the sensing matrix elements can be independently 
selected from a random distribution), only the rows of the sensing matrix can be randomly measured in quantum tomography. 
The sensing matrix was randomly selected from the rows of full measurement matrix (modular Kronecker in Eqn.(7)). Then, 
we get the precise probability distribution by multiplying the sensing matrix and the WF. The problem considered now is to 
reconstruct a state from measured insufficient data. Computer simulation has been conducted to demonstrate the effectiveness 
of the aforementioned scheme. All the recovery processes were implemented in Java and C. For -minimization, we adopt a 
linearized bergman iteration algorithm[ 13]. 

To simulate the sampling process, we assign a number to each row of the full measurement matrix, then computer gener- 
ated a series of random numbers representing the rows we select to form the sensing matrix. For simplicity, the experiment 
error was totally ignored. So in this pseudo pixel representation, WF can be just roughly reconstructed. This is partly due to 
measurement process itself, for only a small fraction of the WF is sensed once, which means the sensing matrix is not dense 
enough. To overcome this obstacle, two approaches might help: finding another sparse representation or using a denser sensing 
matrix(measurement quantum state, for example, coherent state). For example, for some type of mixed state ensemble, WF 
varies along a particular direction when we rearrange the two dimensional WF into a one dimensional object carefully, it is 
sparse in a discrete cosine Fourier representation, which means this basis is its sparse basis. 

In spite of some drawbacks, when we deal with a high dimensional quantum state reconstruction or what we need is just a 
rough knowledge of WF, our reconstruction algorithm would show its advantage, it makes the measurements more economical. 
And any other prior knowledge could also be helpful to recovery process, for example, the symmetry of the WF. By exploring the 
coiTespondence between compressed sensing and tomography, we have shown the possibility of a CS-based WF reconstruction 
method. In addition, most convex optimization algorithms related to CS theory deal with real signal. And in our case, WF is 
real function. But when we need to reconstruct complex signals or images, for example, ambiguity functionCAFl llTill which is 
complex in general, problems appear Here, we mention that G.Zweig proposed an algorithm] 15| processing complex signals, 
which converted the the usual £i -minimization in complex space to a linear programming problem in real space. 

In conclusion, we have described a reconstruction method based on compressed sensing to recover the WF in phase space 
tomography. We regard the Wigner function as an image in phase space. The values of WF on different sites corresponds 
to pixels' gray values. Our numerical simulation has shown good recovery for a sparse discrete WF with usually supposed 
insufficient measurements. Such a method can also be easily extended to reconstruct electromagnetic field's WF with just some 
straightforward revisions. We hope that our scheme can serve as a motivation for other quantum information or phase space 
optics researches. Possible relations to the popular maximum-likelihood reconstruction in quantum tomography still need to be 
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explored. By combination with other reconstruction ideas, we believe a more efficient and precise phase recovery should be 

possible. 
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